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The currently observed accelerated expansion of the Universe suggests that cosmic flow dynamics 
is dominated by some unknown form of dark energy characterized by a large negative pressure. This 
picture comes out when such a new ingredient, beside baryonic and dark matter, is considered as 
a source in the r.h.s. of the fleld equations. Essentially, it should be some form of un-clustered, 
non-zero vacuum energy which, together with (clustered) dark matter, should drive the global 
cosmic dynamics. Among the proposals to explain the experimental situation, the ^'concordance 
modei', addressed as ACDM, gives a reliable snapshot of the today observed Universe according 
to the CMBR, LSS and SNela data, but presents dramatic shortcomings as the "coincidence and 
cosmological constant problems'' which point out its inadequacy to fully trace back the cosmological 
dynamics. On the other hand, alternative theories of gravity, extending in some way General 
Relativity, allow to pursue a different approach giving rise to suitable cosmological models where a 
late-time accelerated expansion can be achieved in several ways. This viewpoint does not require 
to flnd out candidates for dark energy and dark matter at fundamental level (they have not been 
detected up to now), it takes into account only the "observed" ingredients (i.e. gravity, radiation 
and baryonic matter), but the l.h.s. of the Einstein equations has to be modified. Despite of this 
modification, it could be in agreement with the spirit of General Relativity since the only request is 
that the Hilbert-Einstein action should be generalized asking for a gravitational interaction acting, 
in principle, in different ways at different scales. We survey the landscape of f{R) theories of gravity 
in their various formulations, which have been used to model the cosmic acceleration as alternatives 
to dark energy and dark matter. Besides, we take into account the problem of gravitational waves 
in such theories. We discuss some successes of /(ii)-gravity (where f{R) is a generic function of 
Ricci scalar R), theoretical and experimental challenges that they face in order to satisfy minimal 
criteria for viability. 
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I. INTRODUCTION 

Theories of gravity, alternative to Einstein's General Relativity (GR), have been proposed to cure the problems 
of the standard cosmological model and, above all, because they arise in quantizations of gravity. These alternative 
gravitational theories constitute at least an attempt to formulate a semi-classical scheme in which GR and its most 
successful features can be recovered. One of the most fruitful approaches thus far has been that of Extended Theories 
of Gravity (ETGs), which have become a paradigm in the study of the gravitational interaction. ETGs are based on 
corrections and extensions of Einstein's theory. The paradigm consists, essentially, of adding higher order curvature 
invariants and/or minimally or non- minimally coupled scalar fields to the dynamics; these corrections emerge from 
the effective action of quantum gravity [l| . 

Further motivation to modify GR arises from the problem of fully implementing Mach's principle in a theory 
of gravity, which leads one to contemplate a varying gravitational coupling. Mach's principle states that the local 
inertial frame is determined by the average motion of distant astronomical objects This fact would imply that the 
gravitational coupHng here and now is determined by the distant distribution of matter, and it can be scale-dependent 
and related to some scalar field. As a consequence, the concept of "inertia" and the Equivalence Principle have to 
be revised. Brans-Dicke theory (3j| constituted the first consistent and complete theory alternative to Einstein's GR. 
Brans-Dicke theory incorporates a variable gravitational coupHng strength whose dynamics are governed by a scalar 
field non- minimally coupled to the geometry, which implements Mach's principle in the gravitational theory 0, 3, 
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Independent motivation for extending gravity comes from the fact that every unification scheme of the fundamental 
interactions, such as Superstring, Supergravity, or Grand Unified Theories exhibit effective actions containing non- 
minimal coupHngs to the geometry or higher order terms in the curvature invariants. These contributions are one-loop 
or higher loop corrections in the high-curvature regime approaching the full, and still unknown, quantum gravity 
regime jl| • Specifically, this scheme was adopted in the study of quantum field theory on curved spacetime and it was 
found that interactions between quantum scalar fields and background geometry, or gravitational self-interactions, 
yield such corrections to the Einstein-Hilbert Lagrangian |6l]. Moreover, it has been realized that these corrective 
terms are inescapable in the effective action of quantum gravity close to the Planck energy Of course, all these 
approaches do not constitute a full quantum gravity theory, but are needed as working schemes toward it. 

In summary, higher order terms in the invariants of the Riemann tensor, such as , R^^R^^, R^""^^ R^i,ua(3, ROR, 
or R D'^R, and non-minimal coupHng terms between scalar fields and geometry such as </>^i?, have to be added to the 
effective gravitational Lagrangian when quantum corrections are introduced. These terms occur also in the effective 
Lagrangian of string or Kaluza-Klein theories when a mechanism of compactification of extra spatial dimensions is 
used 

Prom a conceptual point of view, there is no a priori reason to restrict the gravitational Lagrangian to a Hnear 
function of the Ricci scalar R minimally coupled with matter |9|]- Furthermore, the idea has been proposed that 
there are no exact laws of physics, in the sense that the effective Lagrangians describing physical interactions could 
be stochastic functions at the microscopic level. This property would imply that local gauge invariances and the 
associated conservation hold only in the low energy limit and the fundamental constants of physics can vary ■ 

Besides fundamental physics motivations, all these theories have been the subject of enormous attention in cosmology 
due to the fact that they naturally exhibit an infiationary behaviour which can overcome the shortcomings of the 
GR-based standard cosmological model. The cosmological scenarios arising from ETGs seem realistic and capable of 
reproducing observations of the the cosmic microwave background (CMB) (ill, [l^. [l3| . It has been shown that, by 
means of conformal transformations, the higher order and non-minimally coupled terms can be related to Einstein 
gravity with one or more scalar fields minimally coupled to9 gravity (l3. fisl. \m. [itI. [l^. 

Higher order terms always appear as contributions of even order in the field equations. For example, the term i?^ 
produces fourth order equations yjl], R OR gives sixth order equations [l§. [20|. RD^R eighth order equations [2l] |. 
and so on. By means of a conformal transformation, any second order derivative term corresponds to a scalar field. ^. 
Fourth-order gravity corresponds to Einstein gravity with one scalar field, sixth-order gravity to Einstein gravity with 
two scalar fields, etc. [18, 22] It is also possible to show that f{R) gravity is equivalent not only to a scalar-tensor 
theory, but also to GR plus an ideal fiuid This feature becomes interesting if multiple infiationary events are 

desired, because an early infiationary stage could select very large scale structures (observed as clusters of galaxies 
today), while a later infiationary epoch could select smaller scale structures (observed as galaxies today) \2^, with 
each infiationary era corresponding to the dynamics of a scalar field. Finally, these extended schemes could naturally 
solve the graceful exit problem bypassing the shortcomings of known infiationary models [l3. [23|. 

In addition to the revision of standard cosmology at early epochs with the concept of infiation, a new approach is 
necessary also at late epochs. ETGs could play a fundamental role also in this context. In fact, the increasing bulk 
of data accumulated in the past few years have nurtured a new cosmological model referred to as the Concordance 
Model. The Hubble diagram of type la Supernovae (hereafter SNela) measured by both the Supernova Cosmology 
Project [25I and the High-z Team [2^ up to redshifts z ~ 1, has been the first piece of evidence that the universe 
is currently undergoing a phase of accelerated expansion. Balloon-born experiments, such as BOOMERanG [IS] and 
MAXIMA {23, have detected the first and second peak in the anisotropy spectrum of the CMB radiation indicating 
that the geometry of the universe is spatially fiat. In conjunction with constraints on the matter density parameter 
coming from galaxy clusters, these data indicate that the universe is dominated by an unclustered fiuid with negative 
pressure, generically dubbed dark energy, which is able to drive the accelerated expansion. This picture has been 
further strengthened by the recent precise measurements of the CMB spectrum obtained by the WMAP experiment 
[29l . [sol . l3l| . and by the extension of the SNela Hubble diagram to redshifts higher than one [s^l- An overwhelming 
fiood of papers has appeared following this observational evidence, presenting a great variety of models trying to 
explain this phenomenon. The simplest explanation is the well known cosmological constant A [33|. Although it is the 
best fit to most of the available astrophysical data [2^, the ACDM model fails in explaining why the inferred value 
of A is so tiny (120 orders of magnitude smaller) in comparison with the typical vacuum energy values predicted by 
particle physics and why its energy density is comparable to the matter density today (the coincidence problem). 

As a tentative solution, many authors have replaced the cosmological constant with a scalar field rolling down its 
potential and giving rise to the model referred to as quintessence [34. l35|. Even when successful in fitting the data. 



^ The dynamics of these scalar fields are governed given by a second order Klein-Gordon-like equation. 
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the quintessence approach to dark energy is still plagued by the coincidence problem since the dark energy and matter 
densities evolve differently and reach comparable values for a very limited portion of the cosmic evolution coinciding 
at the present era. To be more precise, the quintessence dark energy is tracking matter and evolves in the same way 
for a long time. But then, at late times, somehow it has to change its behavior from tracking the dark matter to 
dominating as a cosmological constant. This is the coincidence problem of quintessence. 

Moreover, the origin of this quintessence scalar field is unknown, leaving a great uncertainty on the choice of the 
scalar field potential. The subtle and elusive nature of dark energy has led many authors to look for completely 
different scenarios able to give a quintessential behavior without the need for exotic components. To this end, it is 
worth stressing that the acceleration of the universe only calls for a dominant component with negative pressure, but 
does not tell us anything about the nature and the number of cosmic fluids fiUing the universe. This consideration 
suggests that it could be possible to explain the accelerated expansion by introducing a single cosmic fluid with an 
equation of state causing it to act like dark matter at high densities and dark energy at low densities. An attractive 
feature of these models, usually referred to as Unified Dark Energy (UDE) or Unified Dark Matter (UDM) models, is 
that such an approach naturally solves, al least phenomenologically, the coincidence problem. Interesting examples 
are the generalized Chaplygin gas [sll, the tachyon field [131 and the condensate cosmology jH]. A different class of 
UDE models has been proposed |33] in which a single fluid is considered: its energy density scales with the redshift in 
such a way that the radiation-dominated era, the matter era, and the accelerating phase can be naturally achieved. 
These models are very versatile since they can be interpreted both in the framework of UDE models and as a two- 
fluid scenario with dark matter and scalar fleld dark energy. The main advantage of this approach is that a suitable 
generalized equation of state can be always obtained and observational data can be fltted. 

There is a yet different way to address the problem of the cosmic acceleration. As stressed in [5^, it is possible 
that the observed acceleration is not the manifestation of another ingredient of the cosmic pie, but rather the flrst 
signal of a breakdown of our understanding of the laws of gravitation in the infrared limit. From this point of view, 
it is tempting to modify the Friedmann equations to see whether it is possible to flt the astrophysical data with 
models comprising only standard matter. Interesting examples of this kind are the Cardassian expansion [4l| and 
DGP gravity ji^l. In the same framework it is possible to flnd alternative schemes in which a quintessential behavior 
is obtained by taking into account effective models coming from fundamental physics and giving rise to generalized 
or higher order gravity actions (see [i^l for a comprehensive review). For instance, a cosmological constant 
term may be recovered as a consequence of a non-vanishing torsion fleld, leading to a model consistent with both 
the SNela Hubble diagram and Sunyaev-Zel'dovich data of galaxy clusters [45]. SNela data could also be efflciently 
fltted including higher order curvature invariants in the gravitational Lagrangian [i^ . |47| . These alternative models 
provide naturally a cosmological component with negative pressure whose origin is related to the cosmic geometry, 
thus overcoming the problems Hnked to the physical signiflcance of the scalar fleld. 

The large number of cosmological models which constitute viable candidates to explain the observed accelerated 
expansion is evident from this short overview. On the one hand, this overabundance of models signals the fact that 
only a limited number of cosmological tests are available to discriminate between competing theories and, on the 
other hand, it shows that we are facing an urgent degeneracy problem. It is useful to remark that both the SNela 
Hubble diagram and the angular size-redshift relation of compact radio sources |i48i] are distance-based probes of 
cosmological models, so systematic errors and biases could be iterated. From this point of view, it is interesting 
to search for tests based on time-dependent observables. For example, one can take into account the lookback time 
to distant objects since this quantity can discriminate between different cosmological models. The lookback time is 
observationally estimated as the difference between the present age of the universe and the age of a given object at 
redshift z. Such an estimate is possible if the object is a galaxy observed in more than one photometric band since its 
color is determined by its age as a consequence of stellar evolution. It is thus possible to get an estimate of the galaxy 
age by measuring its magnitude in different bands and then using stellar evolutionary codes to choose the model that 
best reproduces the observed colors. 

Coming to the weak-fleld-limit approximation, which essentially means considering Solar System scales, ETGs are 
expected to reproduce GR which, in any case, is flrmly tested only in this Hmit [4§|. This fact is a matter of debate 
since several relativistic theories do not reproduce exactly the Einsteinian results in the Newtonian approximation but, 
in some sense, generalize them. As flrst noticed by Stelle [H^l, an i?^-theory gives rise to Yukawa-like corrections in the 
Newtonian potential. This feature could have interesting physical consequences; for example, certain authors claim to 
explain the flat rotation curves of galaxies by using such terms [53:] . Others I5j| have shown that a conformal theory 
of gravity is nothing but a fourth-order theory containing such terms in the Newtonian limit. Besides, an apparent, 
anomalous, long-range acceleration in the data analysis of the Pioneer 10/11, Galileo, and Ulysses spacecrafts could 
be framed in a general theoretical scheme by taking into account corrections to the Newtonian potential [55| . 

In general, any relativistic theory of gravitation yields corrections to the Newtonian and post-Newtonian (PPN) 
potentials {e.g., [HHl) which test the theory Furthermore, the newborn gravitational lensing astronomy [s^l is 

generating additional tests of gravity over small, large, and very large scales which soon will provide direct measure- 
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ments for the variation of the Newtonian coupling [53|, the potential of galaxies, clusters of galaxies and several other 
features of self-gravitating systems. Such data, very likely, will be capable of confirming or ruling out the physical 
consistency of GR or of any ETC. In summary, the general features of ETGs are that the Einstein field equations 
are modified in two ways: i) the geometry can be non-minimally coupled to some scalar field, and/or ii) higher than 
second order derivatives of the metric appear. In the first case we deal with scalar-tensor theories of gravity; in the 
second case we have higher order theories. Combinations of non-minimally coupled and higher order terms can emerge 
as contributions to effective Lagrangians; then we have higher order-scalar-tensor theories of gravity. 

From the mathematical point of view, the problem of reducing generalized theories to an Einstein-like form has been 
extensively discussed. Under suitable regularity conditions on the Lagrangian and using a Legendre transformation on 
the metric, higher order theories take the form of GR in which one or more scalar field(s) source of the gravitational 
field (see, e.g., [i,[59, 60, 61]). O n the other hand, as discussed above, the mathematical equivalence between models 
with variable gravitational coupling and Einstein gravity has been studied using suitable conformal transformations 
[63.f63|. A debate on the physical meaning of these conformal transformations seems to be ongoing (j6^ and references 
therein). Several authors claim a physical difference between Jordan frame (higher order theories and/or variable 
gravitational couplings) since there is experimental and observational evidence suggesting that the Jordan frame is 
better suited for matching solutions and data. Others state that the true physical frame is the Einstein one according 
to the energy theorems [61|. However, the discussion is open and no definitive conclusion seems to have been reached. 
The problem becomes more involved at the semiclassical and quantum level, and should be faced from a more general 
point of view — the Palatini approach to gravity could be useful to this goal. 

The Palatini approach to gravitational theories was first introduced and analyzed by Einstein himself (6^ . but was 
named as a consequence of an historical misunderstanding [63, [fill . 

The fundamental idea of the Palatini formalism is to consider the torsion-free connection F^^ entering the definition 
of the Ricci tensor, to be independent of the spacetime metric g^j/. The Palatini formulation of the standard Einstein- 
Hilbert theory turns out to be equivalent to the purely metric theory. This property follows from the fact that the 
field equations for the connection r[^^, considered to be independent of the metric, produce the Levi-Civita connection 
of the metric gf_i„. As a consequence, there is no reason to impose the Palatini variational principle instead of the 
metric variational principle in the Einstein Hilbert theory. However, the situation changes if we consider the ETGs, 
which depend on functions of the curvature invariants (such as f{R) theories) or couple non-minimally to some scalar 
field. In these cases the Palatini and the metric variational principles provide different field equations and the theories 
thus derived differ [61, 69]. The relevance of the Palatini approach for cosmological applications in this framework 
has been recently demonstrated [i^ . H . [70| . 

From the physical point of view, considering the metric g^jy and the connection F^^, as independent fields means 
to decouple the metric structure of spacetime and its geodesic structure (the connection F^^, in general, is not the 
Levi-Civita connection of g^^i/). The causal structure of spacetime is governed by g^jy while the spacetime trajectories 
of particles are governed by F^^. This decoupling enriches the geometric structure of spacetime and generalizes the 
purely metric formalism. This metric-affine structure of spacetime is naturally translated, by means of the Palatini 
field equations, into a bi-metric structure of spacetime. Besides the physical metric g^i,, another metric appears. 
This new metric is related, in the case of f{R) gravity, to the connection. The connection F^^ turns out to be the 
Levi-Civita connection of and provides the geodesic structure of spacetime. 

For non-minimally coupled interactions in the gravitational Lagrangian in scalar-tensor theories, the new metric 
g^i, is related to the non-minimal coupling; can be related to a different geometric and physical aspect of the 
gravitational theory. Thanks to the Palatini formalism, the non-minimal coupling and the scalar field, entering the 
evolution of the gravitational fields, are separated from the metric structure of spacetime. The situation mixes when we 
consider the case of higher order-scalar-tensor theories. Due to these features, the Palatini approach could contribute 
to clarify the physical meaning of conformal transformations [t^I • 

In this review paper, without claiming for completeness, we want to give a survey on the formal and phenomenolog- 
ical aspects of ETGs in metric and Palatini approaches, considering the cosmological and astrophysical applications 
of some ETG models. The layout is the following. The field equations for generic ETGs are derived in SeclUl Specif- 
ically, we discuss metric. Palatini and metric-affine approaches. In Sec. [Till the e quivalence of metric and Palatini 
f{R) gravities with Brans-Dicke theories are discussed. In Sec. |lV]we introduce theoretical and experimental viability 
of /(i?)-gravity. Briefiy, we discuss on the correct cosmological dynamics and on the instabilities for a particular case 
of f{R). After we discuss the precence of ghost fields and the wealk field limit for metric approach. Finally we con- 
sider the growth of cosmological perturbations and the Chauchy problem. Cosmological applications are considered 
in Sec lWVll We show that dark energy and the dark matter can be addressed as "curvature effects", if ETGs (in 
particular f{R) theories) are considered. We work out some cosmological models comparing the solutions with data 
coming from observational surveys. As further result in Sec. IVIII , we show that also the stochastic cosmological 
background of gravitational waves can be "tuned" by ETGs. This fact could open new perspective also in the problems 
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of detection of gravitational waves which should be investigated not only in the standard GR-framework. Discussion 
and conclusions are drawn in Sec lVIIII 



II. THE THREE VERSIONS OF f{R) GRAVITY 

In this survey we focus on f{R) gravity (see [tS] for a more comprehensive discussion and a list of references, and 
[t^ I for short introductions to the subject). In these theories the Einstein-Hilbert action^ 

SEH^^jd^x^gR+S^^^ (1) 

is modified to 

S=^Jd^xV^fiR) + S^"^K (2) 

where f{R) is a non-linear function of its argument and S*^™^ is the matter part of the action. Actually, there are two 
variational principles that one can apply to the Einstein-Hilbert action in order to derive Einstein's equations: the 
standard metric variation and a less standard variation dubbed Palatini variation. In the latter the metric and the 
connection are assumed to be independent variables and one varies the action with respect to both of them, under 
the important assumption that the matter action does not depend on the connection. The choice of the variational 
principle is usually referred to as a formalism, so one can use the terms metric (or second order) formalism and 
Palatini (or first order) formalism. However, even though both variational principles lead to the same field equation 
for an action whose Lagrangian is linear in R, this is no longer true for a more general action. Therefore, it is intuitive 
that there will be two version of /(i?)-gravity, according to which variational principle or formalism is used. Indeed 
this is the case: /(i?)-gravity in the metric formalism is called metric f{R)-gravity and /(i?)-gravity in the Palatini 
formalism is called Palatini f{R)- gravity. 

Finally, there is actually even a third version of /(i?) -gravity: metric-affine f{R)-gravity. This comes about if one 
uses the Palatini variation but abandons the assumption that the matter action is independent of the connection. 
Clearly, metric affine /(i?)-gravity is the most general of these theories and reduces to metric or Palatini /(i?)-gravity 
if further assumptions are made. In this section we will present the actions and field equations of all three versions of 
f{R)- gravity and point out their difference. We will also clarify the physical meaning behind the assumptions that 
discriminate them. 

Then brefly has we show above three versions of /(i?)-gravity have been studied: 

• metric (or second order) formalism; 

• Palatini (or first order) formalism; 

and 

• metric-affine gravity. 

These families of theories are discussed in the following. 

A. Metric f{R) gravity 



In the metric formalism the action is 



Smetr^c = / ^^2^ /(i?) + ^^"^ , (3) 



2 Here R is the Ricci curvature of the metric tensor g^^, which has determinant g, G is Newton's constant, and k = SttG. We mostly 
follow the notations of Ref. [t^. 
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and its variation with respect to g^'^ yields, after some manipulations and modulo surface terms, the field equation 

f'{R)R^. - ^ g^. = V^V,/'(i?) - + , (4) 

with a prime denoting differentiation with respect to R, is the covariant derivative associated with the Levi-Civita 
connection of the metric, and □ = V^V^. Fourth order derivatives of the metric appear in the first two terms on the 
right hand side, justifying the alternative name "fourth order gravity" used for this class of theories. 
By taking the trace of eq. (|4]) one obtains 

SnfiR) + Rf{R) - 2f{R) = K T , (5) 

where T = T"a is the trace of the energy-momentum tensor of matter. This second order differential equation for 
f'{R) is qualitatively different from the trace of the Einstein equation R = —kT which, instead, constitutes an 
algebraic relation between T and the Ricci scalar, displaying the fact that f'{R) is a dynamical (scalar) degree of 
freedom of the theory. This is already an indication that the field equations of f[R) theories will admit a larger 
variety of solutions than Einstein's theory. As an example, we mention here that the .Jebsen-Birkhoff's theorem, 
stating that the Schwarzschild solution is the unique spherically symmetric vacuum solution, no longer holds in metric 
f{R) gravity. Without going into details, let us stress that T = no longer implies that i? = 0, or is even constant. 
Eq. ^ will turn out to be very useful in studying various aspects of f{R) gravity, notably its stability and weak- 
field limit. For the moment, let us use it to make some remarks about maximally symmetric solutions. Recall that 
maximally symmetric solutions lead to a constant Ricci scalar. For R = constant and T^^ — 0, eq. ([5]) reduces to 

f{R)R - 2/(i?) = 0, (6) 

which, for a given /, is an algebraic equation in i?. If i? = is a root of this equation and one takes this root, then 
eq. |[4]) reduces to R^^ = and the maximally symmetric solution is Minkowski spacetime. On the other hand, if the 
root of eq. ([6]) is i? = C, where C is a constant, then eq. ^ reduces to R^i, — g^^C /A and the maximally symmetric 
solution is de Sitter or anti-de Sitter space depending on the sign of C, just as in GR with a cosmological constant. 
Another issue that should be stressed is that of energy conservation. In metric f{K) gravity the matter is minimally 
coupled to the metric. One can, therefore, use the usual arguments based on the invariance of the action under 
diffeomorphisms of the spacetime manifold [coordinate transformations x'^ = .x^ + followed by a pullback, 

with the field S,^ vanishing on the boundary of the spacetime region considered, leave the physics unchanged, see 
[tsI to show that T^^ is divergence-free. The same can be done at the level of the field equations: a "brute force" 
calculation reveals that the left hand side of eq. ^ is divergence-free (generalized Bianchi identity) implying that 
V^T'^'' = 0. 

The field equation ^ can be rewritten as form of Einstein equations with an effective stress-energy tensor to the 
right hand side. Specifically, as 



G 



fiiy 



T, 



fJ.V 



(7) 



where 



rieff) 



f{R) - Rf{R) 



■9^.1' 



y^^.nR) 



g^uDfiR) 



(8) 



is an effective energy-momentum tensor constructed with geometric terms. Since T^fJ^'^ is only a formal energy- 
momentum tensor, it is not expected to satisfy any of the energy conditions deemed reasonable for physical matter, in 
particular the effective energy density cannot be expected to be positive-definite. An effective gravitational coupling 
Geff = G/f'{R) can be defined in a way analogous to scalar-tensor gravity. It is apparent that f'{R) must be positive 
for the graviton to carry positive kinetic energy. 

Motivated by the recent cosmological observations, we adopt the spatially flat Friedmann-Lemaitre-Robertson- 
Walker (FLRW) metric to describe the universe. 



(9) 
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where a is the scale factor. Then, the field equations of metric f{R) cosmology become 



2H + 3ij2 = 



(10) 



f'{R) 

f{R)~Rf'{Ry 



pirn) _^ y.///^^) + 2HRf"{R) + Rf"{R) 



(11) 



where H = d/a is the Hubble parameter and an overdot denotes differentiation with respect to the comoving time 
t. The corresponding phase space is a 2-dimensional curved manifold embedded in a 3-dimensional space and with a 
rather complicated structure |7g |. 



B. Palatini f{R) gravity 

In the Palatini version of ,f{R) gravity, both the metric g^i, and the connection are regarded as independent 
variables. In other words, the connection is not the metric connection of g^jy. While in GR the metric and Palatini 
variations produce the same field equations {i.e., the Einstein equations), for non- linear Lagrangians one obtains two 
different sets of field equations.^ 

Palatini f{R) gravity was proposed as an alternative to dark energy, on the same footing as metric f{R) models. 
The original model advanced for this purpose was based on the specific form f{R) = R — n'^/R [70| . 

The Palatini action is 

SpalaUm = ^ [d^xV^ /(i?) + S^'^^ L., V''™^! , (12) 

2k J L J 

where a distinction needs to be made between two different Ricci tensors contained in the theory. R^j^^, is constructed 
from the metric connection of the (unique) physical metric 5^1,, while -R^j/ is the Ricci tensor of the non-metric 
connection r(^^ and defines the scalar R = g^^'^R^^. The matter part of the action does not depend explicitly from 
the connection F^^, but only from the metric and the matter fields, which we collectively label as i/j^™). 
By varying the Palatini action lfT2|) one obtains the field equation 

/'(i?)i?^, - .9^, ^nT^,, (13) 

in which no second covariant derivative of /' appears, in contrast with eq. An independent variation with respect 
to the connection yields 



(v^ nR)g^'') - V, (v^ nmg"^^) <5;^ = o , (i4) 



where denotes the covariant derivative associated to the (non-metric) connection F^^. 
By tracing eqs. lfT3|) and (fT4| we obtain 



f{R)R^2f{R} = KT (15) 

and 

V^(\/^/'(%^'') =0, (16) 

respectively. Eq. iflGjl is interpreted as stating that V-y is the covariant derivative of the "new" metric tensor 

~g^, = f'iR)g^. (17) 



^ By imposing that the metric and Palatini variations generate the same field equations, Lovelock gravity is selected [77l |. GR is a special 
case of Lovelock theory. 



conformally related to g^^i,. Eq. lfT5|) is an algebraic (or trascendental, according to the functional form of f{R)) 
equation for /'(R), not a differential equation describing its evolution. Therefore, f'{R) is a non-dynamical quantity, 
in contrast to what happens in metric f{R) gravity. The lack of dynamics has consequences which are discussed 
below. It is possible to ehminate the non-metric connection from the field equations by rewriting them as 



2(/ 



l\2 



(18) 



C. Metric-afiine f{R) gravity 

The third family of f{R) theories, metric-affine f{R) gravity fld\, is characterized by the fact that also the matter 
part of the action depends explicitly on the connection F, as described by the action 



Sa f fine — j-. 

2k 



(19) 



V^p is possibly a non-symmetric connection, which would lead to torsion associated with matter and to a reincarnation 
of torsion theories. The latter were introduced in view of elementary particles, rather than cosmology, by coupling 
the spin of elementary particles to the torsion. The study of metric-affine f{R) gravity has not been completed yet, 
in particular its cosmological consequences have not been fully elucidated. It is for this reason that our discussion 
will be limited to metric and Palatini f[R) gravity in what follows. 



III. EQUIVALENCE OF METRIC AND PALATINI f{R) GRAVITIES WITH BRANS-DICKE 

THEORIES 

In the same way that one can make variable redefinitions in classical mechanics in order to bring an equation 
describing a system to a more attractive, or easy to handle, form (and in a very similar way to changing coordinate 
systems), one can also perform field redefinitions in a field theory, in order to rewrite the action or the field equations. 

There is no unique prescription for redefining the fields of a theory. One can introduce auxiliary fields, perform 
renormalizations or conformal transformations, or even simply redefine fields to one's convenience. It is important to 
mention that, at least within a classical perspective such as the one followed here, two theories are considered to be 
dynamically equivalent if, under a suitable redefinition of the gravitational and matter fields, one can make their field 
equations coincide. The same statement can be made at the level of the action. Dynamically equivalent theories give 
exactly the same results when describing a dynamical system which falls within the purview of these theories. There 
are clear advantages in exploring the dynamical equivalence between theories: we can use results already derived for 
one theory in the study of another, equivalent, theory. 

The term ' 'dynamical equivalence^ can be considered misleading in classical gravity. Within a classical perspective, 
a theory is fully described by a set of field equations. When we are referring to gravitation theories, these equations 
describe the dynamics of gravitating systems. Therefore, two dynamically equivalent theories are actually just different 
representations of the same theory (which also makes it clear that all allowed representations can be used on an equal 
footing) . 

The issue of distinguishing between truly different theories and different representations of the same theory (or 
dynamically equivalent theories) is an intricate one. It has serious impHcations and has been the cause of many 
misconceptions in the past, especially when conformal transformations are used in order to redefine the fields {e.g., the 
Jordan and Einstein frames in scalar-tensor theory). In what follows, we review the equivalence between metric and 
Palatini f{R) gravity with specific theories within the Brans-Dicke class with a potential. 

Metric f{R) gravity is equivalent to an w = Brans-Dicke theory'' when f"{R) 01, while Palatini modified 
gravity is equivalent to one with uj = —3/2. The equivalence has been rediscovered several times over the years, often 
in the context of particular theories [79| . 



■* The Brans-Dicke action for general values of the Brans-Dicke parameter u) is Sbd = ^ / "^^^ ~ ^ V^t/iV^t/) — V(</') + S'™' . 
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A. Metric formalism 

It has been noticed quite early that metric quadratic gravity can be cast into the form of a Brans-Dicke theory and 
it did not take long for these results to be extended to more general actions which are functions of the Ricci scalar of 
the metric . Let us present this equivalence in some detail. 

We will work at the level of the action but the same approach can be used to work directly at the level of the field 
equations. We begin with metric f{R) gravity. Let f"[R) be non- vanishing and consider the action ([2]); by using the 
auxihary scalar field = i?, it is easy to see that the action 

S^^^jd'x^g[mR-Vm + S^^^ (20) 

with 

i>{4>) = /'(</>) , V{<f) = 0/'(0) - /(</>) (21) 

is equivalent to the previous one. It is trivial that ((2Q|) reduces io ^ li cj) — R. Vice-versa, the variation of ([20|) with 
respect to g^^ yields 

G^^ = ^ (^V^ V - S/.^nV' - g^,^ + ^ ' ^^^^ 
The variation with respect to 0, instead, gives us 

from which it follows that (j) = R because /" 7^ 0. The scalar field (f) = R is clearly a dynamical quantity which obeys 
the trace equation 

3/"(0)n0 + 3/"'(0) WVa0 + 0/'(0) - 2/(0) = kT (24) 
and is massive. Its mass squared 

.2 _ 1 ffo 



""^ " 3 l^f " ^0 j (25) 

is computed in the analysis of small perturbations of de Sitter space (here a zero subscript denotes quantities evaluated 
at the constant curvature Rq of the de Sitter background). It is convenient to consider, instead of 0, the scalar 
ip = f'{(j)) obeying the evolution equation 

iUil; + 2UN^)~ip'^ ^ kT , (26) 
dip 

where f/(^) = F(0(^)) - /((/-(V')). 

To summarize, metric f{R) gravity contains a scalar degree of freedom and the action 

S^^ [ d^x^gm^ Um + , (27) 



2k _ 

is identified as an w = Brans-Dicke theory. This theory ("massive dilaton gravity") was introduced in the 1970's 
in order to generate a Yukawa term in the Newtonian Hmit js^l, and then abandoned. The assumption /" 7^ is 
interpreted as the requirement of invertibility of the change of variable R ipiR)- 

B. Palatini formalism 

In Palatini modified gravity the equivalence with a Brans-Dicke theory is discovered in a way similar to that of the 
metric formalism. Beginning with the action lfT2|) and defining (j) = R and ip = f'{4>), it is seen that, apart from an 
irrelevant boundary term, the action can be rewritten as 



Spalatini — Tj" / d'^^ \/~^ 

Zk j 



2ip 



S-^™) (28) 
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in terms of the metric g^^i, and its Ricci tensor Here we have used the property that, since ^^j^ = g^iv^ the Ricci 
curvatures of g^i, and g^^, satisfy the relation 

^ = i?. + ^v^VV^^-|n^. (29) 

The action (|28l) is easily identified as a Brans-Dicke theory with Brans-Dicke parameter uu = —3/2. 

IV. THEORETICAL AND EXPERIMENTAL VIABILITY OF f{R) GRAVITY 

In order to be acceptable, f{R) theories should not only reproduce the current acceleration of the universe, but 
they must also satisfy the constraints imposed by Solar System and terrestrial experiments on relativistic gravity, and 
they must obey certain minimal requirements for theoretical viability. More precisely, these families of theories must: 

• possess the correct cosmological dynamics; 

• be free from instabilities and ghosts; 

• attain the correct Newtonian and post-Newtonian limits; 

• originate cosmological perturbations compatible with the observations of the CMB and with large scale structure 
surveys; and 

• possess a well-formulated and well-posed initial value problem. 

If a single one of these criteria is not met the theory should be regarded as unviable. In the following we examine 
how f{R) gravity performs with regard to these criteria. 

A. Correct cosmological dynamics 

According to the tenets of standard cosmology, an acceptable cosmological model must contain an early inflationary 
era (or possibly another mechanism) solving the horizon, flatness, and monopole problems and generating density 
perturbations, followed by a radiation- and then a matter-dominated era. The present accelerated epoch then begins, 
possibly explained by /(7?)_^ gravity. The future universe usually consists of an eternal de Sitter attractor, or ends 
in a Big Rip singularity [911. Smooth transitions between different eras are necessary. The exit from the radiation 
era, in particular, was believed to be impossible in many models [sij, but this proved to be not true. In f act, exit 
from the radiation or any era can be obtained as follows. In the approach dubbed "designer f{R) gravity" in [7J|, the 
desired expansion history of the universe can be obtained by specifying the desired scale factor a{t) and integrating 
an ordinary differential equation for the function f{R) that produces the chosen a{t) In general, the solution 

to this ODE is not unique and can assume a form that appears rather contrived in comparison with simple forms 
adopted in most popular models. 

B. Instabilities 

The choice f{R) — R — fJ-'^/R with ^ Hq ^ 10~^^ eV is again the prototypical example model to discuss 
instabilities. Shortly after it was advanced as an explanation of the cosmic acceleration, this model was found to 
suffer from the pernicious "Dolgov-Kawasaki" instability [s^l- This type of instability was later shown to be common 
to any metric fiR) theory with f"{R) < ([84]) and the extension to even more general gravitational theories has 
been discussed |8^. Let us parametrize the deviations from GR as 

f{R) = R + eif{R) (30) 

with e > a small constant with the dimensions of a mass squared and ip dimensionless. The trace equation for the 
Ricci scalar R becomes 



(31) 
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By expanding around a de Sitter background and writing the metric locally as 

g^Lu = Vtii' + hf,^ , (32) 

and the scalar R as 

R=-kT + Ri, (33) 
with Ri a perturbation, the first order trace equation translates into the dynamical equation for Ri 



(f" if" i^p" \e J iip" 

The expression containing dominates the last term on the left hand side, giving the effective mass squared of Ri 

. ^ . (35) 

Therefore, the theory is stable if f"{R) > and unstable if f"{R) < 0. Strictly speaking, GR is excluded by the 
assumption /" 7^ 0, but the well-known stability of this case can easily be included by writing the stability criterion 
for metric f{R) gravity as /" > 0. 

To go back to the example model of [s^l f{R) = R — fJ.'^ / R, this is unstable because /" < 0. The small scale /i 
determines the time scale for the onset of this instability as ~ 10"^^ s [83], making this an explosive instability. 

A physical interpretation of this stability criterion is the following [86]: the effective gravitational coupHng is 
Geff = G/f'iR) and, if dGeff/dR = — f'G j {f'Y > (corresponding to /" < 0), then Ge// increases with R and a 
large curvature causes gravity to become stronger and stronger, which in turn causes a larger i?, in a positive feedback 
loop. If instead dGeff/dR < 0, then a negative feedback stops the growth of the gravitational coupling. 

What about Palatini f{R) gravity? Since this formaHsm contains only second order field equations and the trace 
equation f'{R)R — 2f{R) — kT is not a diflerential equation but rather a non-dynamical equation, as noted above, 
there is no Dolgov-Kawasaki instability [S^l • 

The discussion of metric f{R) instabilities presented above is based on the local expansion l(32|) and, therefore, is 
limited to short wavelength modes (compared to the curvature radius). However, it can be extended to the longest 
wavelengths in the case of a de Sitter background js^l • This extension requires a more complicated formalism because 
long modes introduce inhomogeneities and are affected by the notorious gauge-dependence problems of cosmological 
perturbations. A covariant and gauge-invariant formaHsm is needed here. One proceeds by assuming that the 
background space is de Sitter and by considering the general action 



S = I d^x y/—g 



(36) 



containing f{R) and scalar-tensor gravity as special cases, and mixtures of them. The field equations originating from 
this action become, in a FLRW background space, 

H' = :^(^^' + ^^ + V-3Hf) , (37) 



3/' V 2 2 



-1 



^ + 3iJ0+^fe-|^ + 2^Uo. (39) 



2^ . , 

de Sitter space is a solution of the field equations provided that the conditions 

- /o + 2V^o = , = 2^0' , (40) 

are satisfied. An analysis of inhomogeneous perturbations of small amplitude and arbitrary wavelengths (88| using 
the covariant and gauge-invariant Bardeen-Ellis-Bruni-Hwang formalism j8§| in Hwang's version [9fl] for alternative 
gravitational theories yields the stability condition in the zero momentum limit 

^^°^'75"-^"' > , (41) 
/o/o 
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This is the stability condition of de Sitter space in metric f{R) gravity with respect to inhomoqeneous perturbations 
and coincides with the corresponding stabihty condition with respect to homogeneous perturbations [8a |. 

The equivalence between metric f{R) gravity and an = Brans-Dicke theory holds also at the level of perturba- 
tions; doubts advanced to this regard have now been resolved. The stability condition of de Sitter space with respect 
to inhomogeneous perturbations m uj = Brans-Dicke theory is given again by eq. l(4T|) , while that for stability with 
respect to homogeneous perturbations is 

(/O)' -VofH ^ g ^42) 

/o 

This inequality is again equivalent to (|4T|) if stability against local perturbations {i.e., /g > 0) is also required. Hence, 
metric f{R) gravity and uj = Brans-Dicke theory are equivalent also with regard to perturbations. 

Beyond the linear approximation, metric f{R) theories have been shown to be susceptible to non-linear instability, 
potentially threatening the possibility of constructing models of relativistic stars in strong f{R) gravity. Inside compact 
objects with spherical symmetry, a singularity could develop if R becomes large [92|. Avoiding this singularity requires 
some degree of fine-tuning. Various authors have contended that this problem can be cured by adding, for example, 
a quadratic term aR^ to the action as first [oil. [93|. This problem needs further study, since it could be the biggest 
challenge left for metric f{R) theories. 



C. Ghost fields 



Ghosts are massive states of negative norm which ruin unitarity and appear frequently in attempts to quan- 
tize Einstein's theory. Fortunately, f{R) gravity theories are free of ghosts. More general ETGs of the form 
/ (i?, R^^R^'^, R^i^-faR^^'^" , •■•), in general, are plagued by the presence of ghosts. A possible exception under certain 
conditions studied in |l94] is provided by theories in which the extra terms are restricted to appear in the Gauss-Bonnet 
combination Q — R^ — 4,Ri,,,jR^'^ + R^^^crR^"'^'^ , as in / = / {R,G). Then, the field equations reduce to second order 
equations without ghosts |95l . [9^ . 



D. The weak-field limit for metric f(R) gravity 



After errors and omissions in the early treatments of the weak-field limit of metric and Palatini modified gravity, a 
satisfactory discussion of the particular model f{R) = R — u'^/R in the metric formaHsm appeared jo^l, followed by 
the generalization to arbitrary forms of the function f{R) |98l. [99l|. 

One studies the PPN parameter 7 which is constrained by light defiection experiments in the Solar System. The 
goal consists of finding the weak-field solution of the field equations and, using this solution, computing the parameter 
7. A static, spherically symmetric, non-compact body which constitutes a perturbation of a background de Sitter 
universe is considered, as described by the line element 

ds^ = - [1 + 2*(r) - H^r^] dt^ + [l-f 2^{r) + TJ^r^] dr"^ + r^dfl'^ (43) 

in Schwarzschild coordinates, with dfl^ being the line element on the unit 2-sphere. ^ and <& are post-Newtonian 
potentials with small amplitudes, i.e., |5'(r)| , |$(r)| << 1, and small (non-cosmological) scales such that Hor « 1 
are considered. The Ricci scalar is expanded around the constant curvature of the background de Sitter space as 
R{r) = Rf) + Ri. The PPN parameter 7 is then given by 7 — $(r)/4'(r) [4§|. The analysis rehes upon three 
assumptions [9^ : 

1. f{R) is analytical at Rq; 

2. mr << 1, where m is the effective mass of the scalar degree of freedom of the theory. In other words, this scalar 
field is assumed to be light and with a range larger than the size of the Solar System (there are no experimental 
constraints on scalars with range m^^ < 0.2 mm). 

3. The matter composing the spherical body has neghgible pressure, P ~ and T = Tb + Ti ~ —p. 

While it is easy to satisfy the first and the last assumptions, the second one is more tricky, as discussed below. The 
trace equation ^ turns into 

V2i?i-m2i?i = ^, (44) 
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regulating the Ricci scalar perturbation, where 

= ^-^=7^ (45) 

is the effective mass squared of the scalar, which reproduces the expression derived in the gauge-invariant stability 
analysis of de Sitter space and in propagator calculations. 
If mr < < 1 , the solution of the linearized field equations is 

^ , , —kM 1 , , 

= - ' (46) 



kM 
12^ 



*(0 = TTr-77 - : (47) 



and the PPN parameter 7 is given by 



-$(r) 1 , , 

7 = -TTTZr = o • (48) 



4'(r) 2 

This value manifestly violates the experimental bound |lQCll | 

|7~ 1| < 2.3- 10"^ (49) 

This violation would mark the demise of metric f{R) gravity were it not for the fact that the second assumption 
necessary to perform this calculation is usually not satisfied. In fact, mr fails to be smaller than unity due to 
the chameleon effect. This effect consists of a dependence of the effective mass m on the spacetime curvature or, 
alternatively, on the matter density of the surroundings. The scalar degree of freedom can have a short range (for 
example, m > 10""^ eV, corresponding to a range A < 0.2 mm) at Solar System densities, escaping the experimental 
constraints, a nd h ave a long range at cosmological densities, which allows it to have an effect on the cosmological 
dynamics (osl. Il0l| . While the chameleon effect may seem a form of fine-tuning, one should bear in mind that f{R) 
gravity is complicated and the effective range does indeed depend on the environment. The chameleon mechanism 
is not arranged, but i s bu ilt into the theory and is well-known and accepted in quintessence models, in which it was 
originally discovered | |102| . It has been studied for many forms of the function f{R) which pass the observational 
tests. For example, the model 

2 f R 



f{R)=R-{l-n)ti'[-^] (50) 



is compatible with the PPN limits if /i ~ 10"^*^ eV~ 10^^^ Hq [101| . To understand how this model can work it is 
sufficient to note that a correction ~ i?" to the Einstein-Hilbert Lagrangian R with n < 1 will eventually dominate 
as i? 0+ . The model l(50|) agrees with the experimental data but could be essentially indistinguishable from a dark 
energy model. Discriminating between dark energy and f{R) models, or between modified gravity scenarios should 
be possible on the basis of the growth history of cosmological perturbations. 



E. Growth of cosmological perturbations 

Since the spatially homogeneous and isotropic FLRW metric solves the field equations of many gravitational theories, 
the expansion history of the universe by itself cannot discriminate between various ETGs. However, the growth of 
structures d epen ds on the theory of gravity considered and has the potential to achieve this goal. A typical study is 
that of Ref. |103| : these authors postulate an expansion history a{t) characteristic of the ACDM model and find that 
vector and tensor modes are not affected by f{R) correctio ns to Einstein gravity, to lowest order, and can be neglected, 
whereas scalar modes do depend on the theory chosen. In |l03i | the stability condition f"{R) > discussed above for 
scalar perturbations is also recovered. It is found there that f{R) corrections lower the large angle anisotropies of the 
cosmic microwave background and produce correlations between cosmic microwave background and galaxy surveys 
which are different from those obtained in dark energy models. A rigorous and mathematically self-consistent approach 
to the proble m of cosm ological perturbations in /(i?)-gravity as been developed using covariant and gauge- invariant 
quantities in [loi . [l05l ? ]). 

The study of structure formation in modified gravity is still uncomplete and, most of the times, is carried out within 
specific f{R) models. Insufficient attention has been paid to the fact that some of these models are already ruled out 
because they contradict the weak-field limit or the stability conditions. A similar situation is found in Palatini models 
which, for this reason, will not be discussed here with regard to their weak-field limit and cosmological perturbations. 
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F. The initial value problem 

A physical theory is required to make predictions and, therefore, it must have a well-posed Cauchy problem. GR 
satisfies this requirement for "reasonable" forms of matter [t^. The well-posedness of th e initial value problem for 
vacuum f{R) gravity was briefiy discussed for special metric models a long time ago [l06l |. Owing to the equivalence 
between f{R) gravity and scalar-tensor gravity when f"{R) ^ 0, the initial value problem of /(i?) gravity is reduced 
to the one for Brans-Dicke gra vity with u; = or —3/2. The Cauchy problem was shown to be well-posed for partic ular 
scalar-tensor theories in |1M . Il07l | but a general analysis has been completed only relatively recently [lOSl . Il09l |. A 
separate treatment, however, was necessary for = 0, —3/2 Brans-Dicke theory. 

We begin by defining the basic concepts employed: a system of 3 -I- 1 equations is said to be well-formulated if it 
can be written as a system of equations of only first order in both temporal and spatial derivatives. Assume that this 
system can be cast in the full first order form 

dtu + A'PV.u^ S{u) , (51) 

where u collectively denotes the fundamental variables hij,Kij, etc. introduced below, M* is called the characteristic 
matrix of the system, and S (u) describes source terms and contains only the fundamental variables but not their 
derivatives. Then, the initial value formulation is well-posed if the system of PDEs is symmetric hyperbolic (i.e., 
the matrices M* are symmetric) and strongly hyperbolic if SiM* has a real set of eigenvalues and a complete set of 
eigenvectors for any 1-form Sj , and obeys some boundedness conditions |llCll |. 

To summarize the results of [lll|, the Cauchy problem for metric f{R) gravity is well-formulated and is well-posed 
in vacuo and with "reasonable" forms of matter (i.e., perfect fluid s, scalar flelds, or the Maxwell fleld). For Palatini 
f{R) gravity, instead, the Cauchy problem is well-formulated [ll2l | but not well-posed in general, due to the presence 
of higher derivatives of the matter fiel ds in the field equations and to the fact that it is impossible to eliminate them 
[111]. However, as it was remarked in |ll3l |. the Cauchy problem for Palatini is still well-posed in vacuo and when the 
trace of the matter energy-momentum tensor vanishes or it is a constant. On the other hand, it is possible to sh ow 
the well-formulation and the well-position as soon as the source of the field equations is perfect-fiuid matter |ll3l | . 

As an alternative, the Brans-Dicke theory equivalent to Palatini f{R) gravity can be mapped into its Einstein 
frame representation. In thi s conformal frame the redefined Brans-Dicke field couples min imal ly to gravity and 
non-minimally to matter |ll4| and the non-dynamical role of this scalar is even more obvious [ll4l |. 

The problems with Palatini f{R) gravity manifest themselves from a. com pletely different angle when one tries to 
match static interior and exterior solutions with spherical symmetry jll6|.^ The field equations are second order 
PDEs for the metric components and, since / is a function of R, which in turn is an algebraic function of T due to 
eq. ifTS]) . the right hand side of eq. ifTS)) contains second derivatives of T. Now, T contains derivatives of the matter 
fields up to first order, hence eq. ifTSj) contains derivatives of the matter fields up to third order. This property is 
very different from the familiar situation of GR and most of its extensions, in which the field equations contain only 
first order derivatives of the matter fields. A consequence of this dependence on lower order derivatives of the matter 
fields is that, in these theories the metric is generated by an integral over the matter sources and discontinuities in 
the matter fields and their derivatives are not accompanied by unphysical discontinuities of the metric. In Palatini 
f{R) gravity, instead, the algebraic dependence of the metric on the matter fields creates unacceptable discontinuities 
in the metric and singularities in the curvature, which were discovered in [ll6l |. Both the failure of the initial value 
problem and the presence of curvature singularities with matter fields can be ascribed to the non-dynamical nature 
of the scalar degree of freedom and to the fact that the latter is related algebraically to T. A possible cure consists 
of modifying the gravitational sector of the Lagrangian in such a way that the order of the field equations is raised. 



V. DARK ENERGY AS CURVATURE 



Let us now show, by some straightforward arguments, how /(i?)-gravity can be related to the problem of dark 
energy. The field equations ^ may be recast in the Einstein-like form 

G^. = R^u - \9^.R = Tlfjf^ + T^Jf'iR) (52) 



^ Other problems of Palatini f{R) gravity were reported and discussed in |ll7l . lll8l |. 
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with T'^'^ff^ given by eq. ([8]) and in which matter couples non-minimally to the geometry through the term l/f'{R). 

As noted above, the appearance of f'{R)-p,i, in T^fJ^^ makes eq. (|52| a fourth order equation (unless f{R) = R, in 

which case the curvature stress - energy tensor T^^/-^^ vanishes identically and l(52l) reduces to the second order Einstein 

equation). As is clear from eq. (|52l) . the curvature stress-energy tensor T^t/'^^ formally plays the role of a source in 
the field equations and its effect is the same as that of an effective fluid of purely geometrical origin. However, one 
can also consider the Palatini approach jfil, [ill, in which the Einstein equations can still be rewritten as effective 
Einstein equations containing a fluid of geometric origin. 

In principle, the scheme outlined above provides all the ingredients needed to tackle the dark side of the universe. 
Depending on the scale considered, the effective curvature fluid can play the role of both dark matter and dark energy. 
From the cosmological point of view, in the standard framework of a spatially flat homogenous and isotropic universe, 
the cosmological dynamics are determined by the energy budget through the Friedmann equations. In particular, the 
cosmic acceleration is achieved when the right hand side of the acceleration equation remains positive. In units in 
which SttG = c = 1 this means 

^^~l{ptot + mot) , (53) 

where the subscript tot denotes the sum of the curvature fluid and the matter contributions to the energy density and 
pressure. The acceleration condition a > for a dust-dominated model is 

Peff + PM +3Peff <0 (54) 



or 



Ptot 

Weff < -7, • (55) 

-^Peff 



Then, the effective quantities 

Peff = |i [f{R) - Rf'iR)] - SHRfiR)^ , (56) 

and 

i?/"(i?) + R liifiR) - Hf"{R) 

Wef f = -IH ■■ (57) 

" [f{R)- Rf'{R)]/2-3HRf"{R) 

play a key role in determining the dynamics of the universe. To gain insight into the dynamics, one can begin by 
neglecting ordinary matter and studying the power-law form f{R) = /o-R" (with n a real number), which represents a 
straightforward generalization of Einstein's GR corresponding the n = 1 limit. This choice yields power-law solutions 
for the scale factor a(t) which provide a good fit to the SNela data and are in good agreement with the estimated 
age of the universe in the range 1.366 < n < 1.376 [? ]. The same kind of analysis can be carried out in the presence 
of ordinary matter, but in this case, numerical solution of the field equations is required. Then, it is still possible 
to confront the Hubble fiow described by such a model with the Hubble diagram of SNela. The fit to the data is 
remarkably good (see Fig.[T]) improving the value and it fixes the best-fit value at n = 3.46 if the baryons contribute 
to the energy density by Qb ~ 0.04, in agreement with the prescriptions if Big Bang nucleosynthesis. The inclusion 
of dark matter does not modify the fit appreciably, supporting the assumption that dark matter is not essential in 
this model. From the evolution of the Hubble parameter in terms of redshift, one can even calculate the age of the 
universe tuniv The best-fit value n = 3.46 provides tuniv ~ 12.41 Gyr. Of course, f{R) = foR" gravity represents 
only a toy model generalization of Einstein's theory. Here we only suggest that several cosmological and astrophysical 
results can be well reproduced in the realm of a power-law extended gravity model. This approach allows fiexibility 
in the value of the exponent n, although it would be preferable to determine a model capable of working at various 
scales. Furthermore, we do not expect to be able to reproduce the entire cosmological phenomenology by means of 
a simple power-law model, which is not sufficiently versatile [8l|. For example, it can be easily demonstrated that 
this model fails when it is analyzed with respec t to its ability of providing the correct evolutionary conditions for 
the perturbation spectra of matter overdensities |l20l | . This point is typically regarded as one of the most important 
arguments suggesting the need for darm matter. If one wants to discard this component, it is crucial to match the 
observational results related to the large-scale structure of the universe with the CMB. These carry the imprints of 
the initial matter spectrum at late times and at early times, respectively. It is important that the quantum spectrum 
of primordial perturbations, which provide the seeds of matter perturbations, can be recovered in the framework of 
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Figure 1: The Hubble diagram of twenty radio galaxies together with the "gold" SNela sample is plotted versus the redshift z, 
as suggested in [119]. The best-fit curve corresponds to the /(i?) gravity model without dark matter. 



i?" gravity. In fact, the model /(i?) oc i? + i?^ can represent a viable model with respect to CMB data and is a good 
candidate for cosmological inflation. To obtain th e ma tter power spectrum suggested by this model, we resort to the 
equation for the matter contrast obtained in Ref. 12d | for fourth order gravity. This equation can be deduced in the 
Newtonian conformal gauge for the perturbed metric 



US 



ds^ = - (1 + 2?A) dt^ + (1 + 2(/)) Sf^i(dx')2 . (58) 

In GR, it is (/) = — ■(/' because there is no anisotropic stress; in general, this relation breaks down in ETGs and the 
non-diagonal components of the field equations yield new relations between the potentials and ^. In f(K) gravity, 
due to the non-vanishing with i ^ j, the (j) — ip relation becomes scale-dependent. Instead of the perturbation 

equation for the matter contrast 5, we provide here its evolution in terms of the growth index s = dlnJ/dlna, a 
quantity measured at z ^ 0.15: 



s (a) — 



- + -E\a) 
a a 



1-2Q 3»„,a-4 
s{a) — - — — r • = = 0, (59) 



where E{a) = H{a)/HQ, R is the dimensionless Ricci scalar, and 

_ 2fnRC^ P 

ho- 



rn 



For n = 1, eq. ([GOl gives the ordinary growth index relation of the Standard Cosmological Model. It is clear from 
eq. (|59| that the latter suggests a dependence of the growth index on the scale which is contained in the corrective 
term Q and that this dependence can be safely neglected when Q ^ 0. In the most general case one can resort to 
the Hmit aH < k < lO^^h Mpcr^ in which eq. (|59l) is a good approximation, and non-linear effects on the matter 
power spectrum can be neglected. 

By studying numerically eq. ((59l) one obtains the evolution of the growth index in term of the scale factor. Assuming, 
for simplicity, the initial condition s(a/s) = 1 at the last scattering surface as in the case of matter domination, the 
results are summarized in Fig. [21 which displays the evolution of the growth index in i?" gravity and in the ACDM 
model. 

In the case of flm — ^bar ~ 0.04, one can observe a strong disagreement between the expected rate of the growth 
index and the behavior induced by power-law fourth order gravity models. This negative result is evident in the 
predicted value of 5(02 = 0.15), which has been observationally estimated by the analysis of the correlation function 
for 220,000 galaxies in the 2dFGRS dataset at the survey effective depth z — 0.15. The observational result suggests 
s = 0.58±0.11 [ml, while our model gives 5(0^ = 0.15) ~ 0.117 (fc = 0.01), 0.117 (fe = 0.001), 0.122 (fc = 0.0002). 
Although this result seems frustrating with respect to the underlying idea of discarding the dark components in the 
cosmological dynamics, it does not give substantial improvement in the case of i?" gravity model plus dark matter. 
In fact, it is possible to show that, even in this case, the growth index prediction is far from being in agreement 
with the ACDM model and again, at the observational scale z = 0.15, there is not enough growth of perturbations 
to match the observed large scale structure. In this case one obtains 5(02 = 0.15) ~ 0.29 (fc — 0.01), 0.29 (fc = 
0.001), 0.31 (fc = 0.0002), i.e., values which are substantially increased with respect to the previous case but still 
very far from the experimental estimate. No significantly different results are obtained if one varies the power n 
(of course, for n ^ 1, one recovers the standard behavior if a cosmological constant is added to the model). These 
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Figure 2: The evolution of the growth index / in terms of the scale factor. The left panel corresponds to modified gravity, in 
the case fim = ^bar ~ 0.04, for the SNela best fit model with n — 3.46. The right panel shows the same evolution in the 
ACDM model. In the case of R" gravity it is shown also the dependence on the scale k. The three cases k = 0.01, 0.001, and 
0.0002 have been eamined, and only the last of these three cases revelas a very small deviation from the leading behavior. 



results seem to suggest that an extended gravity model incorporating a simple power-law of the Ricci scalar, although 
cosmologically relevant at late times, is not a viable description of the cos mic e volution at all scales. Such a scheme 
seems too simple to account for the entire cosmological phenomenology. In [12Q|| a gravity Lagrangian considering an 
exponential correction to the Ricci scalar f{R) = R + Aexp(— i? R) (with A, B constants) produces more interesting 
results and exhibits a grow factor rate in agreement with the observational results at least in the dark matter case. To 
corroborate this point of view, one has to consider that when f{R) is chosen starting from observational data in an 
inverse approach as in [? ], the reconstructed Lagrangian is a non-trivial polynomial in the Ricci scalar. This result 
suggests that the whole cosmological phenomenology can be accounted only by a suitable non-trivial function of the 
Ricci scalar rather than a simple power-law. The results obtained in the study of the matter power spectra for simple 
i?" gravity do not invalidate the general approach. 



VI. DARK MATTER AS CURVATURE 



The results obtained at cosmological scales motivate further analysis of f{R) theories from the phenomenological 
point of view. One wonders whether the curvature fluid which works as dark energy could also play the role of 
effective dark matter, providing an opportunity to reproduce the observed astrophysical phenomenology using only 
visible matter (see for a discussion [74;]). It is well known that, in the low energy limit, higher order gravity implies 
a modifled gravitational potential, which will play a fundamental role in our discussion. By considering a spherical 
mass distribution with mass m and solving the vacuum fleld e quati ons for a Schwarzschild-like metric, one obtains 
the modified gravitational potential of the theory f{R) = foR"' |l22l | 



^(r) 



Gm 



(61) 



where 



12^2 -In -I 



12n3 - 83n2 + hQn+l 



6n^ 



4n + 2 



(62) 



This potential corrects the ordinary Newtonian potential with a power-law term. The correction becomes important 
on scales larger than Vc and the value of this threshold constant depends essentially on the mass of the system. The 
corrected potential l(6T|) reduces to the standard Newtonian potential <& oc 1/r for n — 1 , as follows from the inspection 
of eq. ((62l) . 

The result l(6T|) deserves some comments. As discussed in detail in [l22l |. we have assumed spherical symmetry in 
the the weak-field approximation of the field equations, which leads to a corrected Newtonian potential due to the 
strong non-Hnearity of higher order gravity. Note that Birkhoff's theorem does not hold, in general , in f(R) gravity, 
and that spherically symmetric solutions different from the Schwarzschild one exist in these ETGs |l23| . 

The generalization of eq. (|6T|) to extended sources is achieved by dividing the latter into infinitesimal mass elements 
and integrating the potentials generated by these individual elements. An i nteg ral over the mass density of the system 
is calculated, taking care of possible symmetries of the mass distribution [l22l |. Once the gravitational potential has 
been computed, one can evaluate the rotation curve v'^{r) and compare it with the astronomical data. For extended 
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Figure 3: Best-fit theoretical rotation curve superimposed to the data for the LSB galaxy NGC 4455 (left) and NGC 5023 
(right). To better show the effect of the correction to the Newtonian gravitational potential, we report the total rotation curve 
Vc{r) (solid line), the Newtonian one (short dashed) and the correction term (long dashed). 



systems, one typically must resort to numerical techniques, but the main efTect may be illustrated by the rotation 
curve for the point-Hke situation, that is, 



n, , Gm 



,.0 



(63) 



In comparison with the Newtonian result v1 = Gm/r, the corrected rotation curve is modified by the addition of 
the second term on the right hand side of eq. ((63|) . For < /3 < 1, the corrected rotation curve is higher than 
the Newtonian one. Since measurements of the rotation curves of spiral galaxies signal circular velocities larger than 
predicted by the observed luminous mass and Newtonian potential, the above result suggests the possibility that the 
modified gravitational potential of fourth order gravity may fill the gap between theory and observations without the 
need for additional dark matter. 

The corrected rotation curve vanishes asymptotically as in the Newtonian case, while it is usually claimed that 
observed rotation curves are fiat (i.e., asymptotically constant). Actually, observations do not probe Vc up to infinite 
radii but only show a flat rotation curve (within the uncertainties) up to the last measured point. The possibility that 
Vc goes to zero at infinity is by no means excluded. In order to check observationally this result, we have considered a 
sample of low surface brightness (LSB) galaxies with well measured HI and Ha rotation curves extending far beyond 
the visible edge of the system. LSB galaxies are known to be ideal candidates to test dark matter models because 
of their high gas content, which allows the rotation curves to be well measured and corrected for possible systematic 
errors by comparing 21 cm HI line emission with optical Ha and [Nil] data. Moreover, these galaxies are supposed 
to be dominated by dark matter, so fitting their rotation curves without this elusive component would support ETGs 
as alternatives to dark matter. 

Our sample contains fifteen LSB galaxies with data on the rotation curve, the surface mass density of th e gas 
component and i?-band disk photometry extracted from a larger sample selected by de Blok & Bosma (124 • We 
assume that stars are distributed in an infinitely thin and circularly symmetric disk with surface density S(r) = 
Ti,/oexp(— r/rd), where the central surface luminosity Iq and the disk scale length rd are obtained from fitting to the 
stellar photometry. The gas surface density has been obtained by interpolating the data over the range probed by 
HI measurements and extrapolated outside this range. When fitting the theoretical rotation curve, there are three 
quantities to be determined, namely the stellar mass-to-light (M/L) ratio T*, and the theory parameters {(i,rc). It 
is worth stressing that, while fit results for different galaxies should provide the same value of /?, rc is related to 
one of the integration constants in the field equations. As such, this quantity is not universal and its value must be 
determined on a galaxy by galaxy basis. However, it is expected that galaxies with similar mass distributions have 
similar values of Tc so that the scatter in Tc must refiect the scatter in the circular velocities. In order to match the 
model with the data, we perform a likeHhood analysis for each galaxy, using as fitting parameters /3, logrc (with Tc 
in kpc) and the gas mass fractio n^ fg . As it is evident from the results of the different fits, the experimental data are 
successfully fitted by the model |122l |. In particular, from a purely phenomenological point of view and leaving aside 



This is related to the M/L ratio by T* = [(1 - fg)Mg]/(fgLa), where w Mg = I.AMhi is the gas (HI + He) mass, and = T^L^ 
and -Ljj = 27r/o''^ are the disk total mass and luminosity, respectively. 
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for the moment other viability criteria, from the best fit range (3 = 0.80 ± 0.08, one can conclude that i?" gravity 
with 2.3 < n < 5.3 (best fit value n = 3.2 which overlaps well the above-mentioned range of n fitting SNela Hubble 
diagram) can be a good candidate for solving the missing matter problem in LSB galaxies without dark matter. 

At this point, one wonders whether a link may be found between f?" gravity and the standard approach based 
on dark matter haloes since both theories fit equally well the same data. As a matter of fact, it is possible to 
define an effective dark matter halo by imposing that its rotation curve equals the correction term to the Newtonian 
curve induced by i?" gravity. Mathematically, one can split the total rotation curve derived from i?" gravity as 
^c('') = Ar('') + corri''') t where the second term is the correction. Considering, for simplicity a spherical halo 
embedding a thin exponential disk, we may also write the total rotation curve as v'^{r) — w^.^sfel'') + DAii'^) with 
w^^jgj,(r) the Newtonian disk rotation curve and v'^^DA-ii''') — GM£)M{'r)/r the dark matter one, M^Mir) being its 
mass distribution. Equating the two expressions yields 

Af ( \ M f \ 2P-'r,-0{l-(3)7j^Io{7j)-Vd{rj) 

MdmW = Myir — 3 flTs (64) 

with 77 = r/rd and Sq = TJo, Vdiv) = Io{ii/2)Ko{ii/2) x /i (77/2)^^1 (7^/2). ^ Moreover, 

Io{v,P)^ Mv,v\P}k''-''v'^e~'^'dr/ , (65) 



^0 

where JFq depends only on the geometry of the system and the subscript "wzr" indicates virial quantities. Eq. l(64|) 
defines the mass profile of an effective spherically symmetric dark matter halo whose ordinary rotation curve provides 
the part of the corrected disk rotation curve resulting from the addition of the curvature correction to the gravitational 
potential. Clearly, from a phenomenological point of view there is no way to distinguish this dark halo model from 
i?" gravity. 

Having assumed spherical symmetry for the mass distribution, it is straightforward to compute the mass density 
for the effective dark halo as poA-iir) — {1/ 4:TTr'^)dA'lDM / dr . The most interesting feature of the density profile 
is its asymptotic behavior quantified by the logarithmic slope aoM = dlnpoM/dlnr, which can be computed only 
numerically as a function of 77 for fixed values of (or n). As expected, aoM depends explicitly on f3, while (rc. So 7 fd) 
enter indirectly through rjyir. The asymptotic values at the centre and at infinity (ao and ttoo, respectively) are of 
particular interest, ao almost vanishes and, in the innermost regions, the density is approximately constant. Indeed, 
Qfo = is the value corresponding to models with an isothermal sphere as the inner core. It is well known that 
galactic rotation curves are typically best-fitted by cored dark halo models. Moreover, the outer asymptotic slope lies 
between —3 and —2, values typical of most dark halo models in the literature. In particular, for = 0.80 one finds 
(ao,Q!oo) — (—0.002,-2.41), values which are quite similar to those in the Burkert model, (0,-3). This empirical 
model has been proposed to fit the LSB and dwarf galaxies rotation curves. The values of (ao, aoo) found for the best- 
fit effective dark halo therefore suggest a possible theoretical motivation for Burkert-like models. By construction, the 
properties of the effective dark matter halo are closely related to the disk properties, hence some correlation between 
the dark halo and the disk parameters is expected. In this regard, exploiting the relation between the virial mass and 
the disk parameters, one obtains the relation for the Newtonian virial velocity Vvir = GMvir/rvir 



" 2/^^(1 -/3)G^ MK-W) • ^''^ 



We have checked numerically that eq. l(66|) may be well approximated by Md oc FJJ^ . This relation has the same formal 
structure of the baryonic Tully-Fisher (BTF) relation Mf, oc Vj^i^^ where Mi, is the total (gas plus stars) baryonic mass 
and Vfiat is the circular velocity on the fiat part of the observed rotation curve. In order to test whether the BTF can 
be explained by the effective dark matter halo proposed, we should look for a relation between Vvir and Vfiat- Such 
a relation cannot be derived analytically because the estimate of Vfiat depends on the peculiarities of the observed 
rotation curve, such as how far it extends, and the uncertainties on the outermost points. For given values of the disk 
parameters, we simulated theoretical rotation curves for some values of Vc and measured Vfiat finally choosing the 
fiducial value for Tc that gives a value of Vfiat as close as possible to the measured one. Inserting the relation thus 
found between Vfiat and Vmr into eq. ((66l) and averaging over different simulations, we finally obtain 

logMb^ (2.88±0.04)logy/iat + (4.14±0.09) , (67) 



Here /; and Ki, with I = 1,2 are the Bessel functions of first and second type, respectively. 
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while a direct fit to the observed data gives [l25l | 



logAff, = (2.98 ±0.29) log y/,at + (3.37 ±0.13) . 



(68) 



The slope of the predicted and observed BTF are in good agreement, lending further support to our approach. The 
zero point is markedly different from the predicted one, being significantly larger than the observed one. However, 
both relations fit the data with a similar scatter. A discrepancy in the zero point can be due to our approximate 
treatment of the effective halo which does not take into account the gas component. Neglecting this term, we should 
increase the effective halo mass and hence Kir which affects the relation with Vfiat leading to a higher than observed 
zero point. Indeed, the larger Mg/M^, the more the points deviate from our predicted BTF thus confirming our 
hypothesis. Given this caveat, we can conclude with some confidence that R" gravity offers a theoretical foundation 
even for the empirically found BTF relation. 

Although the results outlined here pertain to the simplistic choice f{R) — foR^ of fourth order gravity, they are 
nevertheless interesting. The incompatibility of this model with the correct matter power spectrum suggests that a 
more complicated Lagrangian is needed to reproduce the entire dark sector phenomenology at all scales, but it has been 
shown that ETGs allow one to approach important issues in cosmological and astrophysical phenomenology. We have 
seen that ETGs can reproduce the SNela Hubble diagram without dark matter and predict the age of the universe. The 
modification of the gravitational potential arising in higher order gravity could constitute a fundamental ingredient in 
interpreting the fiatness of the rotation curves of LSB galaxies. Furthermore, if one considers the model parameters 
selected by the fit of the observational data of LSB rotation curves, it is possible to construct a phenomenological 
analog of the dark matter halo with shape similar to that of the Burkert model. Since the latter has been empirically 
introduced to account for the dark matter distribution in LSB and dwarf galaxies, this result provides a theoretical 
motivation of the Burkert model. 

By investigating the relation between dark halo and disk parameters, a relation has been deduced between Md and 
Vfiat, which reproduces the baryonic Tully-Fisher law. Exploiting the relation between the virial mass and the disk 
parameters, one can obtain a relation for the virial velocity which can be satisfactorily approximated as Md cc T^lr- 
This result is also intriguing because it provides a theoretical interpretation of another phenomenological relation. 
Although not definitive, these phenomenological aspects of f{R) point to a potentially interesting avenue of research 
and support the quest for a unified view of the dark side of the universe. 



As we have seen, a pragmatic point of view could be to "reconstruct" the suitable theory of gravity starting from 
data. The main issues of this "inverse " approach is matching consistently observations at different scales and taking 
into account wide classes of gravitational theories where "ad hoc" hypotheses are avoided. In principle, as discussed 
in the previous section, the most popular dark energy cosmological models can be achieved by considering f{R) 
gravity without considering unknown ingredients. The main issue to achieve such a goal is to have at disposal suitable 
datasets at every redshift. In particular, this philosophy can be taken into account also for the cosmological stochastic 
background of gravitational waves (GW) which, together with CMBR, would carry, if detected, a huge amount of 
information on the early stages of the Universe evolution. In this section we discuss the cosmological background of 
gravitational waves (GWs) in generic f{R) theories. The achievement of detecting massive modes or selecting f{R)- 
signatures in the stochastic background could be the final way to retain or rule out such theories with respect to GR. 
GWs are perturbations ft-^i/ of the metric which transform as 3-tensors. In GR, the equations ruling the propagation 
of GWs in the transverse-traceless gauge are 



where Latin indexes run from 1 to 3. We want to derive the analog of eq. l(69l) for a generic f{R) theory described by 
the action The hnearized theory in vacuo (5^™) — 0) is considered below, so that 



Using a conformal transformation, the scalar degree of freedom /'(i?) of metric f{R) gravity appears as the conformal 
factor in 



VII. MASSIVE SCALAR MODES OF f{R) GRAVITATIONAL WAVES 



(69) 




(70) 



(71) 
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The conformally equivalent Einstein-Hilbert action is 



+ ^ 



where £ ($,$5;^) is the scalar field Lagrangian obtained using the relation 



R ^ e~2* (R ~ 6n$ - 6$.A<I>''') 



(72) 



(73) 



between the Ricci curvatures of the conformally related metrics g^u and g^^- The equation for the gravitational waves 
is now 



dM = , 



where 



Since scalar and tensor modes are decoupled, we have 

H = g^'Sgu = e-'^g'^e'-^Sga = ^ , 



(74) 
(75) 
(76) 



which means that hi is conformally invariant. As a consequence, the plane wave amplitudes hi [t) — h(t)el exp(iA;ma;™), 
where el is the polarization tensor, are the same in both metrics, a fact that is important in the following. 
In a FLRW background, eq. l(74|) becomes 

'SH + 2^^ h + k^a-^h^O (77) 

where k is the wave number and h is the amplitude. The solutions of this equation are linear combinations of Bessel 
functions. Several primordial mechanisms generating GWs are possible. In principle, one could seek for contributions 
due to all known high-energy processes in the early phases of the cosmic history. 

Here we consider the background of GWs generated during inflation, which is strictly related to the dynamics 
of the cosmological model. In particular, one can assume that the main contribution to this background comes 
from the amplification of vacuum fiuctuations at the transition between the infiationary phase and the radiation era. 
However, we can assume that the GWs generated as zero-point fluctuations during inflation undergo adiabatically 
damped (~ 1/a) oscillations until they reach the Hubble radius H~^. This is the particle horizon for the growth of 
perturbations. Any previous fluctuation is smoothed away by the inflationary expansion. The GWs freeze out for 
a/k:S> and re-enter the horizon after reheating. The re-entry in the Friedmann era depends on the spatial scale 
of the GWs. After re-entry, GWs are in principle detectable due the Sachs- Wolfe effect that they induce on the CMB 
temperature anisotropy AT/T at decoupling. If <& is the inflaton fleld, then $ ^ during inflation. By using the 
conformal time 77 defined by drj ^ dt/a, eq. ((77l) becomes 



h" + 2—h' + k^h = Q. 
X 



(78) 



where x = ci^d ^ prime now denotes differentiation with respect to ry. Inside the radius , it is fcry ^ 1. Since 
there are no gravitons in the initial vacuum state, only negative-frequency modes appear and the solution of eq. l(78|) 
is 



h = k^/^ J2hi-^C eyi^ii-ikri) 
all 



(79) 



where C is the amplitude parameter. At the first horizon crossing aH = k, the averaged amplitude Ah — (fc/27r)'^/^ \h\ 
of the perturbation is 



(80) 



When the scale a/k becomes larger than the Hubble radius H^^, the growing mode freezes. It can be shown that 
the upper limit AT/T < is valid since other effects can contribute to the background anisotropy. From this 
consideration, it is clear that the only relevant quantity is the initial amplitude C in eq. l(79|). which is conserved 
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until re-entry. This amplitude depends on the fundamental mechanism that generates the perturbations. Inflation 
produces perturbations as zero-point energy fluctuations, a mechanism which depends on the gravitational interaction 
and (AT/T) further constrains the theory of gravity. Considering a single graviton in the form of a monochromatic 
wave, its zero-point amplitude is obtained from the canonical commutation relation 

[h{t,x),Trh{t,y)]^i5^x-y) (81) 

at flxed time t, where the amplitude h is the fleld and tt/j is the conjugate momentum operator. The Lagrangian for 
the /i-quantity is 

C=lV^r''h;f,h.,, (82) 

in the conformally rescaled FLRW metric g^,^, where the amplitude h is conformally invariant. This Lagrangian leads 
to 

Trh = ^ = e^-^a^h (83) 
dh ^ ' 

Mt,x), /i(y,y)l . (84) 



and eq. l(8T|) becomes 



The flelds h and h can be expanded in terms of creation and annihilation operators. The commutation relations in 
conformal time are 

[hh'*^h*h']^'-^. (85) 

Eqs. ^ and ^ yield C = V^ir^He-'^, where H and $ are calculated at the flrst horizon crossing and, using 
e^* = /'(i?), the relation 

H 

Ah = , 86 

is found to hold for a generic f{R) theory. This result deserves some discussion. Clearly, the GW amplitude produced 
during inflation depends on the theory of gravity which, if different from GR, contains extra degrees of freedom 
which could be probed by the Sachs- Wolfe effect. This effect could be combined with other constraints on the GW 
background if ETGs are probed independently at other scales [l26l . Il27j . 
We are by now familiar with the trace of the fleld equations 

3Df'{R) + Rf'iR) ~ 2f{R) = 0, (87) 

and, using the identiflcations [l28l | 

$^/'(i?) and % ^ V{R)-Rf'{R) (gg) 
the Klein-Gordon equation for the effective scalar fleld <& 

» = ^ (89) 

follows. Lineari zing around a constant curvature background corresponding to <& = $o, assuming that V has a 
minimum at $o |128| . and expanding as in 

V^^a6^\ ^^m2<5$];, (90) 

where the constant m has the dimensions of a mass, yields 
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Figure 4: The evolution of the GW amplitude for a few power-law choices of the scale factor a{t) ~ , the scalar field ~ t™, 
and the function f{R) ~ -R". The horizontal (time) and vertical (amplitude) scales depend on the cosmological background 
providing a signature of the model. 

to first order in h^^i, and (5$. If R^vpa, Rfiv, and R are the linearized quantities corresponding to Rfj.up(j, Rp.v, and i?, 
then the linearized field equations are 



'Oihf — rrP'hf , 



where 



(5$ 

The curvature tensor R^^^pa and eqs. IpTj) are left invariant under the gauge transformations 
By introducing 



(91) 



(92) 



(93) 



(94) 



and considering the gauge vector given by 



(95) 
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the Lorenz gauge 

df'h^, = (96) 
can be chosen. In this gauge the field equations assume the form 

□V = . (97) 

ahf:^m^hf. (98) 
The solutions of eqs. (|97l) and l(98|) are the plane waves 

hf,„ = Afj,^{~p) expiip'^Xa) + c.c. , (99) 

hf — a{'J)') exp{iq°'xa) + c.c. , (100) 

with 

fc" = ([j,'p) uj=p=\'jj\ 

(101) 



q°' = {uj„i, ~p) uj,n = \J rii^ + ■ 

Eqs. I|97p and i(99|) are the wave equation for standard GR and its gravitational wave solutions, respective ly, where as 
eqs. l(98|) and (jlOOp are the wave equation and its solution for the massive scalar mode of /(-R) gravity (cf. [l28lll29l |). 
The dispersion relation for the modes of the massive field hf is non-linear. "Ordinary" (i.e., GR) tensor modes h^y 
propagate at the speed of li ght c , but the dispersion law (the second of eqs. IjlOip ) for the scalar modes /i/ is that of 
a massive field wave packet [123 . Il29l |. The group velocity of a wave packet of /i/ centered on ~p is 



^ = — , (102) 

which is the velocity of a massive particle with mass m and momentum 'p . The second of eqs. (|10ip in conjunction 
with eq. (fT02l) yields 



and a wave packet propagates at constant speed if 



VG = (103) 



• (104) 



The Lorenz gauge is preserved by gauge trasformations with Dsi, = 0; this gauge imposes the transversality condition 
k^Ai_i_y = for the tensor modes, but not for the field h^^ which contains a scalar mode, as seen from eq. l(94|) . or 



h^j^t, h^j^u - ^ri^,y + ri^yhf. (105) 
Were the scalar mode massless, one could impose that 

□e^ = , 

(106) 

thus obtaining a transversal "total" field. However, as is clear from the previous sections, we are dealing with a massive 
scalar mode and transversality is impossible. By applying d'Alembert's operator to the second of eqs. (|106p and using 
eqs. ((97]) and it follows that 

□e^ = rn^hf , (107) 
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in contrast with the first of eqs. I|106p . Similarly, it is shown that a linear relation between the tensorial modes ft-^i^ 
and the massive scalar hf cannot exist. Thus, a gauge in wich /i^i, is purely spatial cannot be chosen, i.e., it is 
impossible to impose hfj,o — 0, see eq. I|105p . However, the traceless gauge condition can be imposed on /i^j,. 



□e^ = , 



(108) 



implying that 



The gauge transformations 



V = ■ (109) 



□e^ = , 
5,,e^ = , 



(110) 



preserve the gauge cJ^/i'"' = 0, h — 0. By choosing ~]) along the z-direction, a gauge can be chosen in which only An, 
A22, and A12 — A21 are different from zero, with the condition h = providing An = —^22- The substitution of 
these equations into eq. (|105p then yields 

V(t, z) = A+{t - z)el+J +A''{t- z)ell'^ + hj{t - vgz)i^^, . (Ill) 

The term A'^{t— z)e'^u (t— z)e^^-' describes the two standard polarizations of tensor gravitational waves familiar 

from GR, while the term hf{t — VGz)rj^,y is the massive scalar field characteristic of f{R) gravity. As expected, the 
scalar /' (i?) generates a third massive polarization for gravitational waves which is absent in GR. 

VIII. CONCLUSIONS 

Let us emphasize once more that we regard f{R) gravity theories not as definitive theories, but rather as toy models 
and proofs of principle that modifying gravity at large scales can explain the observed acceleration of the universe 
without the need to advocate exotic dark energy. This hope has stimulated a very intense activity among theoreticians 
([t^I and references therein). 

To summarize the status of modified gravity, let us note that metric f{R) gravity models exist that pass all the 

- 1 



observational and theoretical constraints (see, e.g., the Starobinsky model [103] f{R) = R + XRo 

The viable models require the chameleon mechanism in order to pass the weak-field limit tests. 

All metric f{R) theories must satisfy the condition f"{R) > to avoid the Dolgov-Kawasaki local instability. This 
is a condition on short-wavelength modes. The stability condition (jUJ is valid for arbitrary wavelengths, but is 
restricted to de Sitter space (which is, anyway, an adiabatic approximation for slowly expanding FLRW spaces). An 
important open problem is whether curvature singularities appear, in general, in relativistic strong field stars. 

As far as the Palatini formalism is concerned, the central idea of this version of modified gravity is to regard the 
torsion-free connection F^^ as a quantity independent of the spacetime metric g^i, . The Palatini formulation of the 
standard Hilbert-Einstein theory is equivalent to the purely metric theory, as a consequence of the fact that the 
field equations for the connection give the Levi-Civita connection of the metric q^u- 

Therefore, there is no reason to 

impose the Palatini variational principle, instead of the metric variational principle, in the Einstein-Hilbert theory. 
However, the situation is difefrent in ETGs containing non-linear functions of the curvature invariants, such as f{R), 
or non-minimally coupled scalars. In these cases, the Palatini and the metric variational principle yield different field 
equations and different physics Wu. 63| . The relevance of the Palatini approach for cosmological applications has been 
amply demonstrated [4^, J0|, O] ■ However, Palatini /(i?) theories could have some problems due to the fact that 
they could contain non-dynamical scalar field and the initial value problem could be ill-posed. In any case, when the 
trace of the matter energy-momentum tensor vanishes identically or it is a constant, and when it can be recast in a 
perfect-fiuid form, the Cauchy problem results well-formulated and well-posed. 

Metric-affine gravity has not been developed in sufficient detail to assess its viability according to all the criteria 
presented here, and its cosmological consequences are essentially unexplored. 
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It seems fair to say that f{R) theories of gravity can help to progress in our understanding of the pecuHarities of 
GR in the wider landscape of relativistic theories of gravity. Furthermore, these theories point out important aspects 
of generaHzations of GR, and, from a phenomenological point of view, constitute viable alternatives to dark energy 
models in explaining the cosmic acceleration, an d to dark matter in reproducing dynamical features as the galactic 
rotation curves or the halo of clusters of galaxies [l3Cll |. Finally, it is possible to "tune" the stochastic background of 
GWs and this occurrence could constitute a further cosmological test capable of conf irming or ruHng out ETGs once 
data from interferometers, Hke VIRGO, LIGO and LISA, will be available (see |l3l| for a discussion on this topic). 
At present, no definite prediction sets f{R) theories apart from dark energy and other models once and for all, but it 
is hoped that progress will me made in this direction. 
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